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Abstract 

Following a recent paper by N. Mandache (Inverse Problems 17 (2001), 
pp. 1435-1444), we establish a general procedure for determining the in- 
stability character of inverse problems. We apply this procedure to many 
elliptic inverse problems concerning the determination of defects of vari- 
ous types by different kinds of boundary measurements and we show that 
these problems are exponentially ill-posed. 

1 Introduction 

Many inverse problems associated to partial differential equations concern the 
problem of determining a parameter of the equation, for example either a co- 
efficient of the equation {coefficient identification) or the geometry (that is the 
boundary) of the region where the phenomenon modelled by the equation oc- 
curs {boundary identification). In order to determine this parameter one needs 
additional information on the solutions to the partial differential equation, usu- 
ally constituted of measurements of the solutions on an accessible (and therefore 
known) part of the boundary of the region in which the phenomenon takes place. 

As an example of a coefficient identification problem, we may think of the 
inverse conductivity problem, whose formulation is due to A. P. Calderon, see [H|. 
In this problem, electrostatic measurements of voltage and current are collected 
on the boundary of a conducting body and by these data one tries to obtain 
information about the conductivity inside the body. 

For what concerns boundary identification problems, we consider the fol- 
lowing examples. First, determination of a defect inside a conducting body by 
electrostatic measurements on the boundary. The defect can be of many dif- 
ferent types: it can be an inclusion, that is a region where the conductivity is 
different from the background conductivity, see for instance and jj^j; it can 
be a crack, that is a fracture, as it has been introduced in ^21; see also 0] and 
|19| : it can be a cavity or a boundary material loss, due to corrosion for example, 
see, for instance, P] and |19j . Then, also the determination of an obstacle by 
acoustic measurements in the far-field can be considered as a problem of this 
kind, see [TU| . 

tSupported by MIUR under grant n. 2002013279. E-mail: dicristoamat.unimi.it 
Supported by MIUR under grant n. 2002013279 and under Progetto Giovani Ricercatori. 
E-mail: rondiSmathsunl .univ . trieste . it 



1 



It has been noted several times that these kinds of inverse problems are ill- 
posed; in fact, even if the amount of data collected is sufficient to guarantee 
uniqueness, the coefficient or the defect, respectively, usually does not depend 
continuously, that is in a stable way, from the measured data. 

For the numerical treatment of inverse problems, the ill-posedness consti- 
tutes a severe difficulty. The second main difficulty is usually due to the fact 
that inverse problems are tipically non-linear, even if the direct problem, in 
the examples above a boundary value problem for an elliptic partial differen- 
tial equation, is a linear one. An accurate knowledge of the character of the 
ill-posedness is an advantage for devising efficient numerical methods. Since the 
problem is ill-posed, in order to recover some kind of stability, we need to apply 
a regularization procedure, that is to restrict the space of admissible unknowns 
(either the coefficients or the defects) by assuming that they satisfy a priori 
conditions involving usually some kind of smoothness assumptions. With this a 
priori information, it is possible to prove that the unknowns depend in a con- 
tinuous way from the measured data. However, an explicit knowledge of such a 
continuous dependence is crucial for several reasons. First, it provides us with a 
quantitative information on how much ill-posed the inverse problem is, thus how 
much difficult it is to solve it numerically; second, a precise knowledge of the 
modulus of continuity of the dependence of the unknowns from the measured 
data indicates the optimal rate of convergence for regularization schemes and 
can be useful also for tuning the regularization parameter, see for instance [llj . 

The determination of the modulus of continuity has to be done in two steps. 
First, we have to establish stability estimates conditioned to some a priori as- 
sumptions on the unknowns; second, we have to show that these stability esti- 
mates are optimal or at least essentially optimal. In order to fulfil the second 
part of this program, we need to construct examples which show that the in- 
verse problem has an instability character of the same order, or at least of the 
same kind, that is of logarithmic or Holder type, for instance, of the stability 
estimates already established. We say that our inverse problem is exponentially 
ill-posed, or severely ill-posed, if such a modulus of continuity is of logarithmic 
type. In other words, exponential instability corresponds to the fact that optimal 
stability estimates are at most of logarithmic type. 

The first of these examples has been constructed in 2 and deals with the 
problem of the determination of a boundary material loss in a planar conductor. 
This example shows that the stability estimates developed in [T^ are essentially 
optimal; since these estimates are of logarithmic type, this kind of problem is 
therefore exponentially ill-posed. An example similar to the one in 2 has been 
constructed for the problem of cavities, still in two dimensions, in 7 . These two 
examples are explicit in the sense that a family of solution showing the instability 
character of the problem is given by explicit formulas, choosing defects whose 
boundaries are highly oscillating. The construction of a family satisfying the 
instability property looked for is not an easy task for other inverse problems. 

Recently, however, N. Mandache has proved in 17 that the inverse conduc- 
tivity problem is also exponentially unstable, showing at the same time that the 
estimates given in P are optimal. The procedure used in ^T] does not depend 
on an explicit construction, it is instead constituted by a purely topological 
argument, which follows from the work of A. N. Kolmogorov and V. M. Ti- 
homirov, ^Hl- We wish to illustrate the argument as follows. Let F : X i-^ Y 
be a function, X and Y being metric spaces. As a model of an inverse problem. 
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X represents the space of unknowns, Y the space of the measured data and F 
is the forward map representing the direct problem. Let us assume that there 
exists xq & X so that for every £ > the baU B{xf),e) contains /(e) disjoint 
balls of radius e/2, /(e) being an integer depending on e. Furthermore, we as- 
sume that for every 5 > there exists an integer g{5) so that F{X) can be 
covered by g{5) balls of radius 5. If we can find ei > and (5(e) so that for every 
e, < e < £i, /(e) > g((5(e)), then we can find xi and X2 in B{xq,£) so that 
dx{x\,X2) > £ and dY{F{xi), F{x2)) < 2S{e). Thus, S{e) provides an indica- 
tion of the instability of the inverse to the map F. Hence, it appears clear that 
establishing this instability character depends on an accurate counting either of 
the maximal amount of disjoint balls with fixed radius that can be found in a 
given ball of the space X or of the minimal amount of balls with fixed radius 
required to cover the image through F oi X. 

This procedure immediately appears to be very general and very well suited 
to be applied in the context of inverse problems. In fact, the space of unknowns 
has, in general, a richer structure with respect to that of the data, since usually 
in inverse problems the forward map F is compact. A first application of this 
procedure to ill-posed problems is developed in |21j . 

Following the topological arguments of jJSj and the procedure described 
in (l7j ■ we have extracted a general method for determining instability, to be 
applicable to many different inverse problems. In Theorem 13.11 below, we have 
stated in a rather abstract framework the outline of the method, in one of its 
possible formulations (for slightly different but analogous formulations we refer, 
for instance, to the discussion of the inverse scattering case, see Subsection l5.4|l . 

Then we have applied our abstract result to many inverse boundary value 
problems of elliptic type. We have shown that all the kinds of boundary iden- 
tification problems briefly described above are exponentially ill-posed. We also 
wish to remark that, as in the explicit examples of [21 and [7], the ill-posedness 
is of exponential type no matter which and how many measurements we take. 
Our examples in fact deal with the ideal case of performing all possible measure- 
ments. This fact is somewhat surprising since in these boundary identification 
problems a much lesser amount of data is required to have unique identifica- 
tion of a defect and also to have stability estimates; usually a finite number 
of measurements is enough. This shows the difficulty of the problem and that 
performing more measurements or different ones does not solve the problem of 
ill-posedness. 

The plan of the paper is as follows. In order to point out to the reader all 
the inverse problems to which we have successfully applied the method, we first 
describe, in Section[21 all the instability results that are contained in the paper. 
Then we proceed with the proofs of these results. The proofs are divided into 
three sections. In Section|3we state and prove an abstract result, Theorem l3.ll 
which provides the general procedure for obtaining the instability examples and 
therefore constitutes the key ingredient and crucial part of the proofs of all 
the results described in Section [21 In fact, the proofs are in general obtained 
as straightforward applications of this abstract theorem. In order to apply the 
abstract theorem, what is essentially needed is to choose a suitable orthonormal 
basis and to check that all the hypotheses of the abstract theorem are satisfied. 
Concerning orthonormal basis, we shall employ eigenfunctions corresponding 
to eigenvalue problems of Stekloff type. We have collected all the information 
we shall need about these orthonormal basis in Section 0] Then, in Section |S1 



3 



the proofs of the instabihty results are concluded. Using the orthonormal basis 
introduced in Section^ we verify that the abstract result applies to the problems 
we consider and we prove their exponential instability. 

In details, in Section |21 first we need to introduce some notations which will 
be used repeatedly in the paper. In particular we define, and investigate the 
structure of, the metric spaces of the unknowns. Then, we list the problems for 
which we have obtained the instability examples, together with the precise for- 
mulation of the instability results. We observe that, for the sake of brevity, we 
usually refer to the bibliography for a more detailed description of the problems 
considered. We begin with the problem of determination of defects of differ- 
ent types by electrostatic boundary measurements. In Subsection 12.11 we treat 
the problem of determination of an inclusion, in Subsection 12. 21 the determina- 
tion of cracks is considered, in Subsection 12 .31 we deal with the inverse problem 
of cavities, in Subsection 12.41 we treat the case of cracks reaching the bound- 
ary of the domain, that is surface cracks, and in Subsection 12.51 we study the 
problem of a boundary material loss. Finally, in Subsection 12.61 we deal with 
inverse scattering problems, in particular with the determination of obstacles 
(either of sound-soft or of sound-hard type) by far-field acoustic measurements. 
In Section 13 the abstract result is stated and proved. In Section 0] we study two 
different eigenvalue problems of Stekloff type and we investigate the asymptotic 
properties of either their eigenvalues or eigenfunctions, in particular this is done 
for three different domains of our interest where the solutions can be computed 
almost explicitly. In Section |31 the conclusions of proofs of all the instability 
results are developed. 

2 Statement of the instability results 

Before stating the main results, we need to introduce some notations about the 
Sobolev spaces we shall use and to describe the spaces of the unknowns. 

For any N > 2, any x — {xi, . . . ,xn) & and any r > 0, we denote 
BNix,r) = {y e : \\y - x\\ < r}. We set S^-'^{x,r) = dBN{x,r). Fur- 
thermore, we set S^-^ = dBN{0, 1), and S^-'^ = {y e S"^"^ : Vn > 0}, and, 
analogously, S'^^^ = {y G 5^^^ : yN < 0}. Finally, we denote B'j^_-^{x,r) — 
{y G Bn{x, r) : y^ ^ xn}- 

We need, furthermore, to introduce the following definition. 

Definition 2.1 Let (F, dy) be a metric space. For a given positive S, Yi, a 
subset of Y, is said to be a S-net for Y if for every y E Y there exists yi G Yi 
so that dY{y,yi) < S. 

Given e positive, Y2 dY is e- discrete if for any two distinct points y2, y2 in 
Y2 we have dr(2/2, 2/2) ^ ^■ 

Notations on Sobolev spaces 

Let il C M^, iV > 2, be a bounded domain and let dil be its boundary. About 
regularity, we assume that there exists a homeomorphism x ■ Bn(0^1) i— > 
such that, for a positive constant C, we have 

^2 2) Wxix) - xiy)\\ < - y\\ for any i,y G Ba,(0, 1), 

||x~^(a;)-X~^(y)ll < C'l|a;-y|| for any a;, y G rJ. 
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Furthermore, wc shall consider two internally disjoint subsets of dfl, Ta and 
Tj, so that U F/ = dil. We assume either that F^ — d^l and Tj — 0, or that 
Ta and F/ are not empty and are assumed to be regular enough, namely there 
exists a homeomorphism x ■ Bn{0, 1) i-^ H. satisfying (|2.1|l . so that, if we still 
denote with x its extension by continuity to Bn(0, 1), then F^ = x{S^~^) and 

rj = x{s!!''). 

We introduce the following Sobolev spaces. Let H^{^l) — {u ^ i^(il) : Vm G 
i^(ri)}, where Vit denotes the gradient of u in the sense of distributions. We 
recall that H^{Q) is a Hilbert space with scalar product {u,v)H^(n) = Jq Vu • 
Vw + uv. With H^/'^iTA) we denote the space of traces of H^{n) functions on 
Ta, which can be endowed in a canonical way with a scalar product induced 
by the one of ffi(O) so that H'^/'^(Ta) is a Hilbert space. By iJ-i/2(r^) we 
shaU denote the dual space to H^/'^{Ta)- We recah that H'^^'^{Ta) C L^CF^) C 
H^^^^{Ta)- We shall also make use of the following spaces. Let oH^^^{Ta) — 
{tP e i/i/2(FA) : /r^ V' = 0}. Its dual is given by the space oH'^^^i^A) = {v e 
H-^/'^{Ta) ■■ iv, 1) = 0}, where (•, •) denotes the duality panmg. 

If F/ is not empty, we set Hq{Q,Ti) and H^^j^^^{n,Ti) as the closed sub- 
spaces of H^{n) constituted by the functions u S H^{n) so that u = in a 
weak sense on Tj and u = constant in a weak sense on F/, respectively. With 
H^^^{Ta,^) and i/^/^gj(F^, fJ) we denote the closed subspaces of H^/'^{Ta) 
constituted by the traces of HQ{il,Tj) and H^^^^^{^l,Tj) functions on F^, re- 
spectively. 

For our purposes, we need to introduce on the Sobolev spaces defined above 
suitable scalar products, which are different but topologically equivalent to the 
canonical ones. We wish to remark that the definitions of these scalar products 
do not take into account the fact that the spaces H"^/^ and H^^^ are dual one 
to each other. 

For any ^, ip £ H^^^{Ta), we set ^ e H^{Q) as the solution to 



(2.2) 



A-0 = in n, 
■0 = V ouTa, 
i = onF,, 



and ip as the solution to the same boundary value problem with -ip replaced by 
ip, and the scalar product we use on H^^^{Ta) is given by 



(2.3) (V',¥')ffi/^(r^) = / V^-V^+ / ^ip. 



We observe that oH^^^{Ta) coincides with the subspace which is orthogonal, 
with respect to this scalar product, to the constant function 1. 

Any r] S H-^^^{Ta) can be decomposed, in a unique way, into the sum of 
fj, an element of oH^^/^(Ta), and a constant function c{r]). Furthermore, to fj 
we can associate fj g H-^{^1) that solves 




(2.4) 

If, in the same way, we associate to G H^^/^{Ta) the functions (f>, c{(j)) and 
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^, then the scalar product on H ^/^(F^) may be defined as 

(2.5) (r?, <^)H-V2(r^) = f Wfi-W^ + c{rii)c{cj)). 

Jn 

We remark that, with respect to this scalar product, oH~^/^{Ta) is the orthog- 
onal subspacc to the constant function 1. 

1/2 

We take Tj not empty. If belongs to H^^^^^iV a, ^^), then there exist (and 

^1/2 

are unique) G (T a, ^) and a constant function c{il}) so that ip — ^p + c{^p). 
Let tp e H^{Q.) solve 



(2.6) 



1 /2 

Then, if we associate to ip G H^^^g^{rA,^) its corresponding decomposition 

1 /2 

given by Lp and c{(p), and its corresponding function (p, on H^^^g^{rA,^) we 
introduce the scalar product 

(2.7) V>)Hlil^ir^,n) = / V^i • + c(V')c((^). 

1/2 

Such a scalar product obviously induces a scalar product on Hq (F^, O), which 

1 /2 

is the closed subspace of ^^^^^^((r^, ft) orthogonal to the constant function 1. 
Spaces of smooth perturbations of a given set 

We shall consider the following examples. Let us fix integers N > 2 and m > 1 
and positive constants s and p. Let us also fix a; G and r > 0. 

To any real function / defined on B^_-^(a;, r), where S^_^(x, r) = {y G 
BN{x,r) : UN = xn}, we associate its graph, that is graph(/) — {y € 
: UN = f{yi,...,yN-i,XN), {yi,...,yN-i,XN) e B'^_-^{x,r)}, and, as- 
suming / > xn, its subgraph, that is subgraph(/) = {y e : xn < yN < 
f{yi,---,yN-i,XN), {yi,...,yN-i,XN) € B^_i(x,r)}. 

With the notation Xjni3e{B'j~f_i{x, r j) we indicate the space {graph(/) : / G 
Ca iBN-i{x,r)), ||./||c'"(s;,_j(x,r)) < /? and a;Ar < / < Xjv + e} and with 
ymi3e{B'f^_i{x,r)) we indicate the space obtained by taking the subgraphs of 
all the functions belonging to the same class as before. We consider the spaces 
Xmi3e{B'j^_i{x,r)) and Ymi3e{B'^_i{x,r)) as metric spaces with the Hausdorff 
distance. 

To any strictly positive function g defined on S^~^{x,r) = dBN{x,r), we 
denote its radial graph as graph^„^(g) ^ {y £ M.^ : y = x + g{u)) ■ {^^^^) , w € 
S^~^{x,r)} and its radial subgraph as subgrapl\.„^(5) = {y € M.^ : y = 
x + p- (^) , < p < g{uj), Lo e S'^-'jx, r)}. 

Then, with the notation Xmfie{S^ ^{x,r)) we denote the space given by 
{graph^ad(ff) : 9 G C"'{S^-^{x,r)), ||5llc'"(S"-i(x,r)) < P and r < g < 
r + e} and with Ym0e{S^~^{x,r)) we denote the space of radial subgraphs 
of all the functions belonging to the same class used before. Also the spaces 
Xmpe{S^~^{x,r)) and Ym0E{S^~^{x,r)) are metric spaces endowed with the 
Hausdorff distance. 
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It is an easy remark the fact that Xm^j(_B^_-^(a;, r)) and X„i/3e(S'^~^(x, r)) 
are contained in the closed baU, with respect to the HausdorfF distance be- 
tween closed sets, of radius e centred at B'i^_^{x,r) and S'^~^{x,r), respec- 
tively. Analogously, Ympe{B']^-i{x,r)) and Ymj3e{S^~^{x,r)) are contained in 
the closed ball, again with respect to the Hausdorff distance, of radius e and 
centre B'^_^[x, r) and Bjs[(x, r), respectively. Maybe more interesting and signif- 
icant is the fact that the elements of Ymf}e{S'^~^{x,r)) are all compact subsets 
which are star-shaped with respect to a common point x e R^. The deter- 
mination of star-shaped sets is usually considered to be more stable than the 
determination of other kinds of sets. Nevertheless many of our examples show 
that even with a star-shapedness assumption the instability is still of exponential 
type. 

We would like to study properties of e-discrctc sets of X„i/3e(i?^_j^(a;, r)), 
Y^0,{B'j,^^{x,r)) and X,r^pe{S''-\x,r)), Y^p.iS^'-^x.r)). We have the fol- 
lowing proposition. 

Proposition 2.2 Let us fix integers N >2 and m > I and positive constants (3 
and r. We also fix x e R^. Fixed e > 0, let be equal to one of the following 
four metric spaces: X„r(}e{B'j^_j^{x,r)), YmpsiB'j^-iix^r)), X.mi3e{S^~^{x,r)) 
orY^p.iS'^-^x.r)). 

Then, there exists a positive constant Sq, depending on N , m, (3 and r only, 
so that for any e, < e < e^, we can find satisfying the following properties. 
We have that the set is contained in X^; Z^ is s-discrete, with respect to the 
Hausdorff distance; and, finally, Z^ has at least exp(2~^£Q^ i)/'"g-(A''-i)/m'j 
elements. 

Proof. The proof can be obtained, with slight modifications, along the lines of 
the proof of Lemma 2 in 1 7 . □ 

2.1 Inverse inclusion problem 

Let us assume that the domain 51 — -Bjv(0, 1), iV > 2, is occupied by a con- 
ducting body. Let us assume that an inclusion D is present inside the otherwise 
homogeneous conductor; that is, there exist two different positive constants a 
and b and a set D which is compactly contained in (that is D is a compact 
subset of fi) so that the conductivity inside D is constantly equal to a and the 
conductivity outside D, that is in fl\D, is constantly equal to b. For the sake 
of simplicity, we normalize b so that 6=1 and we take a to be positive and 
different from 1. 

The electrostatic potential u inside is a solution to the following partial 
differential equation 

(2.8) div((l + (a-l)xD)Vu) = in ft, 

where xd denotes the characteristic function of the domain D. 

Furthermore, u satisfies a boundary condition on (957 = S^~^ which depends 
on whether we prescribe the voltage V' G H^/'^{dQ) on the boundary or we assign 
the current density 77 e QH^^/^{dVL) on the boundary. Namely, in the first case 
the boundary condition is given by 

(2.9) u = V on dn; 
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in the second case by 



(2.10) = ?7 on / u = 0; 

O'' Jon 

where we have added a normahzation condition. 

We have existence and uniqueness of a (weak) solution for both the boundary 
value problems ^^-^1^ and (|TS|l - lfOHl . 

The inverse problem we consider is the one of recovering the shape and the 
location of an unknown inclusion D, by performing current and voltage mea- 
surements at the boundary, that is either by prescribing voltages and measuring 
the corresponding current densities or viceversa. 

In the literature, a lot of attention has been devoted to the determination of 
-D by a single measurement; in this case the problem has been often referred to 
as the inverse conductivity problem with one measurement. A global uniqueness 
result is still missing, see 15] and its references for a more detailed discussion on 
this topic. We remark that if all possible measurements are performed, then the 
inclusion can be uniquely determined, see However, up to our knowledge, 
even if all measurements are considered, no explicit stability estimate for this 
problem has been established. 

We produce an example showing that, even if we make many measurements, 
actually all possible measurements, the optimal stability for this inverse problem 
is at the best of logarithmic type. We treat the case when voltages are prescribed 
and currents are measured and the case in which current densities are assigned 
and voltages are measured, as well. 

We fix two positive integers, m and N, N > 2, and two positive con- 
stants, (3 and a, a 1. We consider the metric space {X,d) where X = 
^m/3(i/4)('S'^~^(0, 1/2)) and d is the Hausdorff distance. If D belongs to X, 
then we can define the following two operators. 

The operator A{D) : H^/^{dn) ^ H-^/^{dn) is defined as 

where i^, e H^/^{dn), u is the solution to H2.8() - (|2.9() and <^ is any H^{^) 
function whose trace on dU, is equal to ip. Since the operator K{D) associates 
the Dirichlet datum to the corresponding Neumann datum, it is usually called 
the Dirichlet-to-Neumann map. 

Viceversa, the operator M{D) : QH^^^^{dQ) oH^/'^{dQ) is given by 

U{D)7j = u\an 

where rj e QH^^^^{dQ) and u is the solution to (|2.8() - (|2.10() . For the same 
reasons, the map M{D) is called the Neumann-to- Dirichlet map. 

It is easy to show that for any D G X, the maps A(r') and A/'(I?) are linear 
and bounded operators between a Hilbert space and its dual. In the sequel, their 
norms will be always assumed to be the canonical ones as bounded operators 
between Hilbert spaces. We state the instability result. 

Proposition 2.3 We fix integers N > 2 and m > I and a positive con- 
stant [3. We also fix Q < a ^ \. Let {X,d) he a metric space, with X = 
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^^/3(i/4)('S'^ "'^(0, 1/2)) and d being the Hausdorff distance. Then we can find 
positive constants Si and C , which depend on N , m, j3 and a only, so that for 
any e, < e < ei, there exist Di and D2 belonging to X such that 

(2.11) d{D,,BN{0, 1/2)) < e, for any j = 1, 2; d{Di,D2)>e; 
and 

(2.12) WHDi) - A{D2)\\ < Ccxp(-£-(^-i)/(2mAr)-). 



(2.13) \\J\f{Di) - J\f{D2)\\ < Cexp(-£-(^-i)/(2™^)). 

The proof of this proposition is postponed to Subsection l5.1l 

Experimental measurements 

In Proposition 12.31 the inverse inclusion problem is stated to be exponentially 
ill-posed even if we perform all possible measurements of current and voltage 
type at the boundary. It is not surprising, therefore, that the inverse inclusion 
problem is exponentially ill-posed also with respect to measurements which can 
be actually obtained from the experiments. We shall refer to this kind of mea- 
surements as the experimental measurements. The model which we shall follow 
is the one developed in 22 , which we briefly describe, referring to the original 
paper for more details. 

The model is the following. On the boundary of the conductor fi, we attach 
L electrodes. The contact regions between the electrodes and the conductor are 
subsets oidft and will be denoted by e;, / = 1, . . . ,L. We assume that the subsets 
ei, I = 1, . . . , L, are open, connected, with a smooth boundary and so that their 
closures are pairwise disjoint. We remark that we identify any electrode with 
its contact region. A current is sent to the body through the electrodes and 
the corresponding voltages are measured on the same electrodes. For each I, 
I — 1, . . . ,L, the current applied to the electrode e/ will be denoted by // and 
the voltage measured on the electrode will be denoted by Vi . The column vector 
I whose components are / — 1, . . . , L, is a current pattern if the condition 
^f^i Ii — is satisfied. The corresponding voltage pattern, that is the column 
vector V whose components are V;, / = 1, . . . , L, is determined up to an additive 
constant and we always choose to normalize it in such a way that X^^^i ~ 
The voltage pattern depends on the current pattern in a linear way, through an 
Lx L symmetric matrix R which is called the resistance matrix, that is y = RI. 

The following model can be used to determine the resistance matrix R. We 
assume that at each electrode ei, I — 1, . . . , L, a surface impedance is present 
and we denote it with zi . Let us assume that there exists Z > so that for each 
I, I — 1, . . . , L, zi > Z . Let D be as before an inclusion in fl. The conductivity 
in D is a, where a is a positive constant different from 1, and the conductivity 
outside D is equal to 1. If we apply the current pattern / on the electrodes, then 
the voltage u inside the body satisfies the following boundary value problem 

{div((l + {a- l)xi3)Vw) = in Q, 

u + zijT ^ Ui on ei, 1^1,..., L, 

|7 = ondn\{jl,eu 

L.^^Ii forany? = l,...,L, 
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where Ui, I — 1, . . . , L, are constants to be determined. We call U the column 
vector whose components are given by ?7;, Z = 1, . . . , L. 

For any Z = 1, . . . , L, VJ, a component of the voltage pattern V , is given 
hyVi^ /e, u, thus, by dH, 

Vi = \ei\Ui-ziIu 

where \ei\ denotes the surface measure of e;. 

By 1^21 Theorem 3.3], we infer that there exists a unique couple {u,U), 
u being in H^ijl) and U being a column vector with L components so that 
Tlii=i \^iWi — zili = 0, such that (|2.14() is satisfied. Thus the current pattern 
/ uniquely determines the voltage pattern V, if this is normalized in such a 
way that X^z^i = 0- Furthermore, it has been proved in [32] that the relation 
between / and V is linear, thus the resistance matrix R{D) is well defined. 
Finally, it has been shown that R{D) is actually symmetric. We remark that 
we shall assume, without loss of generality, that R{D)[1] — 0, where [1] denotes 
the column vector whose components are all equal to 1. Also, we recall that the 
norm of R{D) will always be the norm of linear operators from into itself. 

The following instability result will be proved in Subsection 15. II 

Proposition 2.4 We fix integers N > 2 and m > 1 and a positive constant 
[3. We also fix < a 1 and Z > 0. Let {X,d) be a metric space, with X = 
^m/3(i/4)(*S'^~^(0, 1/2)) and d being the Hausdorff distance. Let us assume that 
L > 2 electrodes ei, I ~ 1, . . . , L, and their surface impedances zi, I — 1, . . . , L, 
are fixed and satisfy the previously stated assumptions. Then we can find positive 
constants Si and C , which depend on N , m, (3, a, Z and the electrodes only, so 
that for any e, < £ < ei, there exist Di and D2 belonging to X such that 

diD,,BN{0, 1/2)) < e, for any j ^ 1,2; d{Di,D2)>e; 
^ ' \\R{D,) - R{D2)\\ < (7exp(-e-(^-i)/(2™^^)). 



2.2 Inverse crack problem 

Let O — Bn{0, 1), > 2, be the region occupied by a homogeneous conducting 
body. Let us assume that inside the conductor there is a crack a, that is a 
closed set inside $7 so that is connected and, locally, a can be represented 
by the graph of a smooth function. We can consider two different types of 
cracks, perfectly insulating and perfectly conducting, and we can prescribe on 
the (exterior) boundary of either the voltage or the current density. Thus, the 
electrostatic potential u in f2 satisfies either 



(2.16) 



Au = in n\a, 
on da; 



if a is perfectly insulating, or, when a is assumed to be perfectly conducting, 

(2.17) 1 A^ = inf^V, 

I u = constant on a. 

We remark that, in (|2.16(l . on da means on either sides of a. On the boundary 
the potential satisfies either 

du 

(2.18) u = ^ondn; (— |ao,l)=0; 
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where tp £ H^/'^{dfl) is the prescribed vohage at the boundary, or, if we pre- 
scribe the current density on the boundary to be 77 G oH^^^^{dQ), 

(2.19) = 77 on 917; / u ^ 0; 

ov Jan 

we wish to remark that normahzation conditions have been added to the bound- 
ary conditions. 

We have that all the direct problems ^J^-^J^, ^J^-^J^, (|TT7|l - (|0^ 
and llTT7jl - l|ni^ admit a unique (weak) solution. 

The inverse crack problem consists of recovering the shape and location of 
an unknown crack a by performing electrostatic measurements at the boundary. 

In this subsection we shall state the instability character of such an inverse 
problem, in all the possible cases, that is when we consider either insulating or 
conducting cracks, and when either we prescribe voltages and measure corre- 
sponding currents or we prescribe currents and measure corresponding voltages. 

For a detailed analysis of uniqueness and stability of this problem we refer to 
|19j. for the two-dimensional case, and to U, for the three-dimensional case, and 
to their bibliographies. We wish to remark that, for what concerns uniqueness 
and stability results, these have been obtained with a finite number of boundary 
measurements, usually with two suitably chosen measurements. Our instability 
example shows the optimality of the stability estimates previously obtained and 
that the stability can not be improved by taking different or more measurements. 

The framework of our example is as follows. Let N > 2 and m, positive 
integers, and /?, a positive constant, be fixed. Let X = ^m/3(i/4)(i?7v-i('^' 1/2)) 
with the Hausdorff distance. To any a G X, we can associate the following four 
operators. 

Let Ai(ct) : H^''^{dn) ^ H-^/^{dn) be given by 

Ou f 

(Ai(cr)7A,(p) = (— = / Vit-V^, 

where ^p, ip £ H^^^{dVl), u solves H2.16|l - H2.18|l and (p is any H^(Vl\a) function 
whose trace on 90 coincides with Lp. 

Let Ni[(7) : oH-^/^{dn) ^ oH^/^{dn) be given by 

Ni{(j)i-i = 

where G oi?"^^^(90) and u solves ^^TT^ - ^^TV^ . 
Let A2 (cr) : H^/^{dQ) ^ H-^''^{dVl) be given by 

du f 
Jn 

where ip, if € H^^'^{dVL), u solves H2.17ll - H2.18|l and ip is any ^/"^^^^((O, cr) function 
whose trace on 90 coincides with (p. 

Let N2{<7) ■■ oH-^/^{dn) ^ QH^/^{dn) be given by 

A/'2(cr)77 = u\an, 

where 77 e oH-^^'^{dn) and u solves (|?T7ll - (|?ll^ . 

Let us remark that for any a £ X, each Ai{a) and Ni{a), i ~ 1,2, is a 
bounded linear operator between a Hilbert space and its dual, it is self-adjoint 
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and its norm is always assumed to be the canonical one of bounded operators 
between these two Hilbert spaces. Keeping in mind these notations and this 
remark, we are able to state our instability result. 

Proposition 2.5 Let us fix integers N > 2 and m > 1 and a positive constant 
j3. Let {X,d) be the metric space where X = X„^(i/4')(_B^_-^(0, 1/2)) and d is 
the Hausdorff distance. Let us fix T e {Ai, A/i, A2, A/'2}. Then there exists a 
positive El, depending on N , m and (3 only, so that for any e, Q < e < Si, there 
exist two cracks ai, belonging to X satisfying 

d{aj, B'j^_^{0, 1/2)) < e, for any j = 1,2; d{cri,a2) > e; 
^ ' llr(fTi) -r(a2)|| < 2exp(-e-(^-i)/(2™JV)), 

For the proof of this proposition we refer to Subsection 15.21 

2.3 Inverse cavity problem 

The inverse cavity problem can be treated if we substitute, in the previous 
subsection, the set of cracks inside V,, X = Xjni3(i/4){B'j,^_i{0,l/2)), with the 
set of cavities inside given hy X = 5^m/3(i/4){5'^~^(0, 1/2)). With almost no 
modification in the proof, a result completely analogous to the one described in 
Proposition 12.51 can be obtained. So also the inverse cavity problem shows an 
exponential instability character. 

We recall that, concerning the inverse cavity problem, stability estimates of 
logarithmic type have been obtained for the two dimensional case in and 
for the higher dimensional case in [3]. For the planar case, an explicit example 
developed in [7] shows the exponential instability character of the inverse cavity 
problem and, consequently, that the stability estimates therein contained are 
essentially optimal. Our results here confirm this fact and extend it to the higher 
dimensional case, thus providing the essential optimality of the estimates proved 
inEl. 

2.4 Inverse surface crack problem 

Let n = Bn {0,1)\{x e BAr(0, 1) : a;Ar_i > and xn = 0}, TV > 2, and let 
F/ = {a; e B]^{0,1) : x^-i > and xn = 0}. Inside il, we consider the 
geodesic distance between two points as the infinium of the lengths of smooth 
paths contained in Q connecting the two points. If we consider the boundary 
of Q with respect to this distance, we notice that this boundary contains two 
overlapping copies of Tj, one obtained by approaching Tj with points a; in f2 
such that Xn > and the other obtained by approaching it with points a; G f2 
so that a;Ar < 0. The set F^ is obtained from F/ by taking the closure, in the 
topology of dil induced by the geodesic distance defined above, of 9r2\F/. We 
remark that F^ coincides, from a set point of view, with S^~^, but each point 
belonging to the intersection of F/ and S^^^ should be counted with multiplicity 
two, as for points of F/. 

With iji/2(r^) 

we denote the space of traces of H^{n) functions on F^ 
and with H~^^^{Ta) we shall denote its dual. On these two spaces, we consider 
scalar products which are defined exactly as we have done before for regular do- 
mains, in and respectively. We notice that H^/^^Ta) C L^{S'^~'^) C 
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H-^^^{Ta)- Finally, we notice that the spaces oH^^^{Ta) and oH-'^^^{Ta) are 
the orthogonal subspaces, respectively in 7?^/^(r^) and i7~^/^(r^), to the con- 
stant function 1 and are dual one to each other. 

We observe that the spaces HQ{il,Tj) and i/^Q„sj(il, F/) are given by the 
spaces of H^{Bn{0,1)) functions which are, respectively, identically zero or 
constant in a weak sense on Tj. The spaces of traces on F^ of functions belonging 
to HQ{Q,ri) and i?^Q„^((ri, F/), respectively, are again denoted by Fq''^ (Fa, fl) 

1 /2 

and i?^Q„5((FA, On these two last spaces, a scalar product is defined in the 
same fashion as we have done for regular domains in (|2.7|) . 

If iV = 2, let ao = {x e ^2(0, 1) : xi > -1/2 and X2 = 0}. If TV > 3, let 
/ e C^{B'j,_^{0, 1/4)) so that -1/4 < / < 0. Let ao = F/ U {y G B'j^_^{0, 1) : 
/((yi, . . . , yN-2, 0)) < yjv-i < 0, {yi, . . . , yN-2, 0) € B'j,_^{Q, 1/4)}. By defini- 
tion if = 2, and by a suitable choice of / if > 3, we can always assume 
that B']^_^{xq, 1/16) is contained in o-q, where xq = (0, • . . , 0, —1/8, 0). 

Then we fix a positive integer m and a positive constant (3 and we define X 
as the set 

(2.21) X = {a-(ao\B;,_i(io,l/16))Ua': a' e X„^(i/4)(S^_i(io, 1/16))}. 

We remark that each cr G X is a connected closed set inside i3jv(0, 1) so that 
Ti da and cr\F7 C Bjv(0,4/5). 

If we assume that i?Ar(0, 1) is occupied by a homogeneous conductor, we can 
think any cr G X as a surface crack inside Bn{0, 1). We can distinguish between 
two different kinds of surface cracks, namely insulating and conducting. 

Let us assume that a G X is an insulating surface crack and that we prescribe 
on Ta the current density to be equal to 77 G o-ff^^^^(FA)- Then the electrostatic 
potential u in Bn{0, 1) satisfies 

{Aii==0 in BAr(0, l)\cr, 
^ ^ on cither sides of a, 
$U = ^ on Ta, 
Iva " = 0, 

where we have also added a normalization condition. We have that m is a weak 
solution to (|2.22|l if and only if u G H^{Bn{0, l)\cr), /p^ m = 0, and 

/ Vm • Vw — rj{w |r^), for any w G H^{Bn{Q, l)\cr). 

v'Bjv(0,1)\ct 

Clearly such a function u exists and is unique. We have that u|ryi belongs to 
oH^^^Ita) and that the operator Ar3(cr) : pg-V^^rA ) ^ oH'^/'^(Ta) so that, 
forany?7Goi?"'/'(FA),AA3(a)7; = u|r^, u solution to l|2.22|l . is linear, bounded 
and self-adjoint. 

When, otherwise, cr G X is a conducting surface crack in Bff{0, 1) and we 

1 /2 

prescribe the voltage on F^ to be G H^ljj^^f{rA,^), then the potential u in 
Bn{0, 1) solves 

( Au = mBN{0,l)\a, 

(2.23) 1 " "^(^^ °^ ^' 

I w = 7/1 on Ta, 
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1/2 ~ 

where c{'tp) is a constant so that = ip — c{ip) G Hq' (r^,f2). Let -0 be any 
-//q (i?Ar(0, l)\cr, (t) function so that ip\rA — i'- Then u solves in a weak sense 
if and only if w - c[ip) - ip & H^{Bn{0, l)\cr) and 

/ Vw • Vw = 0, for any w G (BAr(0, 1)\ct). 

JBjv(0,1)\o- 

By standard elliptic equations methods we infer that u, solution to l|2.23|l . exists 
and it is unique. To such a solution we can associate ^{fa ^ (-^c/nst (r^ , 

1 /2 

as follows. For any (p £ H^^^g^(^A,^), let the constant c{(p) be so that ip — 
(f — c{lp) e i/Q^^(r^, il), and let (p be any Hq{Bj^{0, 1)\(t, cr) function so that 
<^|r^ = 0- Then, 



The operator A^ia) : iji/i,(r^, 1]) (i?i/l,(r^, 1]))' so that, for any V € 
i?j[/^s((r^, ri), A3((7)'(/; = fijIrAJ solution to (|2.2;^() . is linear, bounded and 
self-adjoint. 

The inverse surface crack problem consists of the determination of an un- 
known surface crack from suitable information on the operator As or A3, re- 
spectively. The operators A/3 and A3 correspond to electrostatic boundary mea- 
surements. Many papers have treated this problem when a finite number of 
measurements is performed, that is when either As (77) is measured for a finite 
number of different 77 or A3 (■!/)) is measured for a finite number of different tp. We 
refer to 19 and its bibliography for a detailed description of the problem in the 
planar case. In 19' uniqueness and stability estimates of logarithmic type are 
established. In 5 the uniqueness issue for the three-dimensional case is treated. 
In the next proposition we show that also this inverse problem is exponentially 
unstable, thus proving the essential optimality of the stability estimates of |19j . 

Proposition 2.6 We fix integers N > 2 and m > 1 and a positive constant (3. 
Let X he the set of closed sets described in (|2.21(l and let (X, d) he a metric space 
with the Hausdorff distance. Let ctq he defined as before. Let us fixT d {A3, A3}. 
Then we can find e\ > 0, that depends on N , m and (3 only, so that for any e, 
< e < ei, there exist two surface cracks ai, a2 belonging to X so that 



d{(7j:(7o) <e, for any j = 1,2; d(cri, 0-2) > e; 
||r((Ti)-T(a2)|| < 2exp(-e-(^-i)/(2™^)). 



We prove this result in Subsection [5 

2.5 Inverse boundary material loss problem 

Let O = {x e Bn{Q.I) : xn > }, and let Va ^ £ 9B jv(0, 1) : xn > 
0} ^ S^^^ and Ti ^ {x e 5^(0, 1) : xat = 0} = B;,_^(0, 1). Fixed a positive 
integer m and a positive constant /?, let 

(2.25) X = {a = TiUa' : a' G r,„^(i/4)(S^v-i(0, 1/2))}. 

Then every cr g AT is a closed subset contained in fl so that ((T\r/) C Bn{0, 4/5). 
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Am = 


in Q\(T. 




on r^, 




on r^, 


It A = 





We assume that ft is the region occupied by a homogeneous conductor and 
cr G X is a boundary material loss, which might be due to a corrosion phe- 
nomenon, for instance. We assume that Ta is an accessible part of the boundary 
of our conductor, whereas T„ — d{n\a)\rA, that is the other part of the bound- 
ary where the material loss occurs, is not. Also in this case we distinguish two 
kinds of boundary material losses, insulating and conducting. In the first case, 
no current passes through F^, the part of boundary of f2\cr which is contained 
in a. In the second case, the voltage is constant on a. More precisely, we have 
that if a is an insulating boundary material loss and if we prescribe the current 
density on Ta to be equal to 77 e oH^^/^(rA), then the electrostatic potential 
u inside fl\a is the unique solution to 



(2.26) 



where the last line is a normalization condition. Otherwise, if a is conducting, 
then the electrostatic potential u in 17 is given by 

{Am = in n\a, 
u = c{ip) on cr, 
u — ip on Ta, 

1 /2 

where tp £ H^Q„s((r^, f2) is the prescribed voltage on Ta and c(V') is a constant 

1/2 

so that '0 — c{tp) belongs to Hq (Ta, 

In the insulating case, for every a £ X, we define Afi{a) : o-ff^^^^(r^) 
qH^/^{Ta) so that for any 77 e QH~^/^(rA), then Afi{a)r] — wjr^, m being the 
unique solution to H2.26|l . We have that N4,{a) is a linear, bounded and self- 
adjoint operator. 

1/2 

In the conducting case, if cr e X, let us define Ki{a) : H'J^^^^^{T a,^) ^ 
{Hlil,,{TA,m as follows. For any ^, ^ G Hlil,^{TA,n), 

dxi f 
(A4(a)7A,^) = ( — |r^,(^) = / Vu-Vif, 

where u solves H2.27|l and (p is any Hl^^g^{^l, a) so that pIfa = ^^^o A4{a) is 
a linear, bounded and self-adjoint operator, for any a € X. 

The inverse problem consists of the determination of the shape and the loca- 
tion of an unknown boundary material loss a from electrostatic measurements 
performed on the accessible part of the boundary, that is Ta- The case of a 
single electrostatic measurement is particularly interesting and uniqueness and 
stability estimates have been obtained for this kind of problem, see JHI ; and its 
references, for the two-dimensional case and |3] and also ^ for the higher dimen- 
sional case. The stability estimates obtained in and |3] are of logarithmic 
type and they are essentially optimal. In two dimensions, this has been shown 
through an explicit example provided in [2]. In the next proposition, proven in 
Subsection 15.31 we confirm that logarithmic stability is essentially optimal in 
any dimension, no matter how many and which measurements we perform. 
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Proposition 2.7 Let N > 2 andm > I be integers and [3 be a positive constant. 
Let X be defined as in (|2.25l) . endowed with the HausdorjJ distance d. Let us fix 
T S {A/Ci, A4}. Then there exists a constant Si > 0, that depends on N, m and 
f3 only, so that for any e, < s < Si, there exist two boundary material losses 
<7i, <J2 belonging to X so that 



d{<J],B'^_^{Q,l)) <£, for any j = 1,2; d{ai,a2)>e; 
||T(ai) -T(a2)|| < 2 exp(-e-(^-i)/(2"^)). 



2.6 Inverse scattering problem 

We turn our attention to inverse scattering problems, in particular to tlie deter- 
mination of obstacles inside a medium by acoustic far-field measurements. For 
a detailed description of this kind of inverse problems we refer to jJUj . 

Let us fix integers N E {2, 3} and m > 1 and two positive constants (3 and 

a. 

Let X — Ymi3(i/2){Sl'^~^{0, 1)). We assume that D G X is an obstacle in an 
otherwise homogeneous medium. 

The incident field is determined by a time-harmonic acoustic plane wave 
and is given by u^{x;Lu;a) = e'"^^", where x € K.^, ^/a is the wave number, 
and CO G S^~^ is the direction of propagation. The direct scattering problem 
consists of finding the total field u{x;w,a), x G R^\D, which is the sum of the 
incident field u^(x;uj;a) and of the scattered field u'^{x;uj;a), which is due to 
the presence of the obstacle D. The total field u satisfies 



(2.29) 



Au + au = mR^\D, 

u{x; uj; a) = e'v^^''^ + u^(x; uj; a) for any x G R^\D, 



with a boundary condition on dD which depends on the nature of the obstacle. 
Namely, if the obstacle is sound-soft, then 

(2.30) u = on dD, 
if the obstacle is sound-hard, then 

(2.31) ^ = on dD. 

ou 

Furthermore, the scattered field satisfies the so-called Sommerfeld radiation con- 
dition 

(2.32) hm r(^-i)/2 i^u'\ = 0, 

r^oo \ or J 

where r — \\x\\ and the limit holds uniformly for all directions x = a;/|ja;|j G 
gN-i^ The Sommerfeld radiation condition characterizes outgoing waves and 
implies that the asymptotic behaviour of the scattered field is given by 

(2.33) u^(x;a;;a) = pp^jyyj |u^(i;w;a) + 

as 11 a; 1 1 goes to oo, uniformly in all directions x = a;/||a;|| G S^^^ . The function 
is called the far- field pattern related to the solution to H2.29|l - H2.30|l - H2.32|l . 
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or to l|2.29|) - (|2.31|) - (|2.32|) respectively, for the direction of propagation lo and the 
wave number ^/a. 

Therefore, if G X is sound-soft, then we denote with As{D) : S^^^ x 
S^^^ X (0,00) i-^ C its Jar-field map, that is, for any i, a; G S^^^ and any 
a > 0, 



where is the far- field pattern related to the solution to H2.29|l - H2.30|l - H2.32|l . 

In an analogous way, assuming I? G X to be sound-hard, we can associate 
to D its far-field map, given by Ah{D) : S^^^ x S^'^ x (0,oo) 1-^ C, so that, 
for any ix,Lu,a) G 5^-^ x S'^-^ x (0,oo), 



where is now the far-field pattern related to the solution to 12.29() - (I2.31() - 

The inverse scattering problem consists of the determination of an unknown 
obstacle D from some information about its far-field map. More precisely, we 
assume that we a priori know whether the unknown obstacle is sound-soft or 
sound- hard. Then, from a, usually partial, knowledge of As{D) or Ah{D), re- 
spectively, we try to determine the obstacle D. We remark that information 
about the far-field map can be collected by performing suitable far-field acous- 
tic measurements. 

In the next proposition we show that also this inverse problem is exponen- 
tially unstable. We recall that stability estimates for the determination of a 
sound-soft obstacle have been obtained by V. Isakov, see |14l I15| . 

Proposition 2.8 Let N G {2,3}. Let m and j be positive integers and (3 he a 
positive constant. Let a and a be two constants so that < a <a. We denote 



and we fix any j real numbers oi, . . . , so that Oi G In for any i, i = I, . . . , j . 
Let {X,d) be the metric space with X = Yjni3(i/2){S^~^{0A)) '^'^'^ ^ ^he Haus- 
dorff distance. Then we can find a constant ei > 0, that depends on N, m, j, 
(3 and Ln only, so that for any e, < e < e\, there exist obstacles Di, D2, -D3 
and D4 all belonging to X such that 



(2.34) 



As{D){x,uj,a) = w^(x;w;a). 



(2.35) 



Ah{D){x, uj, a) = ul^{x; uj; a), 



{ 



[a, a] i/iV = 2; 
(0,a] z/7V = 3; 



(2.36) 



d{Dj,BN{0,l)) < e, for any j = 1,...,4; 
d{Di,D2)>e; d{D3,D^)>e; 




sup ||(A(£'i)- A(I>2))(-,-,a)|!L2(S"-ixs«-i) < 2exp(-e 



aG {fli ,. . . ,aj } 




sup ||(y^,j(i:'3) - A(£'4))(-,-,a)|lL2(siv-ixs"-i) < 2exp(-£ "„J). 



aG{ai ,...,aj } 
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We observe that the result is sHghtly different depending whether N = 2 
or N — 3. If N = 2, the real numbers Ui, i — 1, . . . , j, satisfy ai > a and the 
constant ei depends on a as well. If = 3, we do not need a lower bound for 
the numbers a^, apart from the fact that they are positive, and the constant 
£i does not depend on a. We shall point out during the proof, which will be 
developed in Subsection 15.41 where the hypothesis ai > a becomes necessary. 



3 The abstract theorem 

Let {X, d) be a metric space and let H he a separable Hilbert space with scalar 
product (•,•). As usual we denote with H' its dual and for any v' € H' and any 
V G H we denote by {v' , v) the duahty pairing between v' and v. With C{H, H') 
we denote the space of bounded linear operators between H and H' with the 
usual operators norm. We shall also fix 7 : i?\{0} 1-^ [0, +00] such that 

(3.1) 7(Aw) = 7(w) for any v G H\{0} and any A e M\{0}. 

Let us remark that the function 7 may attain both the values and +00 and 
can be thought of as a suitable kind of Rayleigh quotient. 

Let F be a function from X to C{H, H'), that is, for any x ^ X, F{x) will 
denote a linear and bounded operator between H and H' . We also fix a reference 
operator Fq € C{H^ H') and a reference point xq in X. We wish to point out 
that no connection is required between and F^, in particular Fq might be 
different from F{xo)- For any e > 0, lei X^ — {x ^ X : d{x,xo) < e}. 

Recalling the notations introduced in Definition 12.11 we can formulate the 
following exponential instability result related to the map F. 

Theorem 3.1 Let us assume that the following conditions are satisfied. 

i) There exist positive constants Sq, Ci and ai such that for any e, < e < Eq, 
we can find an e-discrete set contained in X^ with at least exp(Cie^"^) 
elements. 

ii) There exist three positive constants p, C2 and oli and an orthonormal basis 
in H , {ffe}^^, such that the following conditions hold. 

For any k € N, we have that j{vk) < 00, and for any n G N, 

(3.2) #{fceN: "fivk) < n} < C2{1 + n)P 

where # denotes the number of elements. 

For any x X and any (/c, ^) G N x N we have 

(3.3) \{{F{x)-Fo)vk,vi) \ < C2cxp{-a2max{-,{vk)n{vi)})- 



Then there exists a positive constant e\, depending on Eq, C\, C2, ai, 012 
and p only, so that for every e, < e < Si, we can find X\ and X2 satisfying 

xi,X2eXe; d{xi,X2)>e; 

^ ' mx^) - F{x2)\\ciH,H') < 2cxp(-£-"i/2(f+l)). 
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A crucial point in the proof of Theorem 13.11 is constituted by the following 
lemma in which we construct (5-nets in the image through F oi X with a control 
in terms of S of the number of their elements. 

Lemma 3.2 Under assumptions^ of Theorem 13. IL there exists a positive con- 
stant C3, depending on p, C2 and ai only, such that for every 6, < 5 < l/e, 
we can find a S-net Yg for F{X) with at most exp(C3(— logfJ)^^"*"^) elements. 

Proof. Let G be an element of C{H, H') and let {wfc}^J?[ be the orthonormal 
basis in H defined in assumption Im}) of Theorem 13. II For any pair (fc, Z) e N x N 
we associate to G the real number ak,i — {Gvk,vi). Let HGjly — supf, ; |afc^;|(2 + 
max{7(t>fe), 7(t)j)})P+^ and let Y be the normed space 

Y = {GeC{H,H') : \\G\\y < 00} 

with norm || • \\y- First, we notice that, for any x G X, F{x) — Fq is contained 
in Y . This is an immediate consequence of (13. 3|) : in fact \{{F{x) — Fo)vkTVi)\ < 
C2 exp(— a2 niax{7(ufc), 7(wi)}) and hence 

\\F{x) - FoWy < sup G2 exp(-a2(n ~ 1))(2 + n)P+^ < 00. 

n 

1 /2 

Second, if we set C4 - C2 {J2„{1 + n)~^) ' , for any G G F we have 

(3.5) WGWciHM') < C4\\G\\y. 

This follows from the following computation 

- yl^k,l\"'k,l\ {2+ma.x{'y{vk),fivi)})2p+^ J - 

^ (Lfc,; (2+max{7(t,fc),7(fi)})^''+V H^"^- 

/ \ 1/2 

Let us show that C4 is a bound for {^^ ^ (2+max{7(,;O.7(t„)})^''+0 ' 
positive integer n, the number of pairs (fc, V) so that n — \< max{7(7jfc), 7(w/)} < 
n is bounded by C|(l + n)^^, therefore 

\ < "S^ + r7^2p \ ^ pi 

(2 + max{7(z;fc),7(z;i)})2p+2 -Z^^2l + ; (2 + (n - l))2p+2 ^4- 

Let us now fix 5, < (5 < l/e. Let n be the smallest positive integer so that 
for any real number i > n it holds C2exp(-a2(t - 1))(2 + <)p+i < (5/(2C4). 
There exists a positive constant G5 , depending on p, C2 and 0L2 only, such that 
n < C5 log (5-1. 

Let S' = ^"^^^l"^'' S and let *5 = S'Z n [-C2, G2]. We remark that -^s is a 
finite subset of R and we have that #^'5 < Cq/5' , where Gg can be chosen as 
2G2 + 1. 

Let us define the following subset of C{H,H'). Let fa = {G G C{H,H') : 
ak,i G ^5 if T^SLx{l{vk) , < n and Uk^i — otherwise}. 
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Let us count the number of elements of Ys . If we set 

s = #{{k,l): max{7(i;fc),7(w;)} < n} 
then we have that #1^ ~ {^'^sY and hence 

#Ys < (2^6^4(2 + n)P+^~^)Cl{l+n?^ < 

< (2^6^4(2 + Cs log<5^i)?'+i<5-i)'='2'(i+^= log 5-^)^- 

and then by straightforward computations we can find a positive constant C3, 
depending on p, C2 and 02 only, so that #Yi < exp(C3(— log(5)^P+^). 

Then we need to show that for every x E X there exists G E Ys so that 
\\Fix)-Fo~G\\ciH,H')<5/2. 

We fix X e X and we set bk,i — {{F{x) — Fo)vk,vi); we recall that F{X) — 
Fq C Y and, by (|3.5|) . it is enough to determine G E Ys so that ||-F(a;) ~ Fq — 
G\\y < 5/{2Ci). We observe that, by 1231, buj G [-C2, C2] for every (fc,/). We 
construct such a G as follows. We set ak,i = {Gvk,vi). For any {k,l) so that 
max{7(ufc), 7(vi)} < n, we prescribe a^^i to be the element of ^'5 that is closest 
to bk,i- If, otherwise, (fc, Z) is so that max{7(ufc), j{vi)} > n, then we set ak,i — 0. 
We have that G belongs to Yg by construction and we can evaluate ||-F(a;) — 
Fq — G\\y as follows. For any {k, I) so that max{7(wfc), j{vi)} < n, \ak,i — bkj\ < 

S', that is \akj - 6fc,i|(2 + max{7K), 7(^0})"+! < (2 + n)P+^^^^%^S < 
5/{2Ci). If {k,l) is such that t = max{7(?;fc), 7(vi)} > h, then ak,i = and 
\bk,i\{2 + t)P+^ < C2exp(-a2t)(2 + t)P+i <5/{2G4) by the definition of n. 

Having established this property, it is easy to find a subset Yg of F{X) with 
the same number of elements as Yg which is a (5-net for F{X) and hence the 
proof is concluded. □ 
Proof of Theorem 13.11 The proof is obtained by combining assumption 0) 
of Theorem 13.11 with Lemma 13.21 as follows. Let e satisfy < e < eo and 5 
satisfy < (5 < 1/e. Let C X^ be, as in assumption [|)| of Theorem l3.1l an e- 
discrete set with at least exp(Ci£~"i) elements. The same procedure employed 
in Lemma [3.21 allows us to find Yg C F{Xir) which is a 5-net for F{X^) with 
at most exp(C3(— log(5)-^^'+^) elements. If > #1^, then there exist xi and 
X2 € X^ so that d{xi^ X2) > s and ||i^(xi) — F{x2)\\c{h,H') < 25. 

In order to have that > #1^5, it suffices to have that exp(Cie~"^) > 
exp{C3(- \og5fP+^). Let us define 5{e) = exp(-e-"i/2(P+i)). Then there exists 
a constant £1, < £1 < Eq, depending on Sq, Ci, C2, ai, a2 and p only, 
such that for every e, < e < ei, we have (5(e) < 1/e and exp(Ci£^"i) > 
exp(C3(— log(5(e))^''+^) and so the result follows. □ 

Remark 3.3 We wish to remark that it is easy to see that the order of insta- 
bility can be improved up to exp(— £""1/*^^^+^+"=*^), for any > 0. However, 
in this case, the constant ei depends on 03, too. For the sake of simplicity, we 
have stated the theorem when is chosen to be equal to 1. 

4 SteklofF eigenvalue problems 

In this section we collect some results which will be repeatedly used later, when 
we shall apply the abstract theorem to find instability examples for inverse 
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problems. Most of the results described in this section arc obtained by standard 
methods, thus, for the sake of brevity, we do not enter into any detail and we 
limit ourselves to fix the notation and to state the results which will be needed 
later, referring to the literature when necessary. 

Let C M^, iV > 2, be a bounded domain and let Ta and F/ be two 
internally disjoint subsets of dfl, so that Ta UTj = dfl. About the regularity 
and the properties of fl, Ta and F/, we shall consider the same assumptions 
used at the beginning of Section 13 at page0] 

The following eigenvalue problems of Stekloff type will be discussed; first 

{Aw = in fl, 
^ = Xu on Fa, 
S; = onF,, 

and then, assuming Tj not empty, 
(EP2) 

We state the following propositions concerning the eigenvalues and eigen- 




functions of ljEPl(l and ljEP2(l respectively. 

Proposition 4.1 Under the assumptions on Q, Ta and F/ previously made, we 



have that the eigenvalues of ljEPl|l . counted with their multiplicity, are given by 
an increasing sequence 

= Ai < Az < A3 < . . . < Afe < . . . 

so that limfc_^oo At = 00. For any n ^ N, we set Ni{n) = i^{k G N : Aj, < n}. 
Then the asymptotic behaviour of the eigenvalues is as follows. There exists a 
constant Ci depending on f2, T a and F/ only so that 

(4.1) Ni{n) <Cin^^^, foranyneN. 

Moreover, there exists a corresponding sequence of eigenf unctions, {uk}keN, 
that is Uk (z H^{n)\{0} and the couple (Afc,Ufc) solves (|EP1|I for any fc £ N, so 
that the following three conditions holds 

{w/clr^lfceN is an orthonormal basis of L'^{Ta)] 
< . ^ \ya \ is an orthonormal basis of H^^^{Ta)', 

{l|r^}U {\/AfcMfe|r^}fc>2 is an orthonormal basis of H~^^'^{Ta)', 

where we have considered the spaces H^^^{Ta) and H-^/^{Ta) with the scalar 
products defined in (|2.3() and H2.5|l respectively. We remark that ui is a constant 
function not identically equal to zero. 

Proposition 4.2 Under the assumptions on fJ, T a and F/ previously made. 



and assuming that F/ is not empty, then the eigenvalues of l)EP2p . counted with 
their multiplicity, constitute an increasing sequence 

< ^1 < ^2 < ■ • ■ < Mfc < • ■ ■ 
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so that linife^oo /^fc = oo. For any n €z N, we set as before N2{n) = =ff{k G N : 
fik < n}. Then the eigenvalues satisfy the following asymptotic behaviour. There 
exists a constant C2 depending on fl, Ta and Tj only so that 

(4.2) N2{n) < C2n'^-\ foranyneN. 

Furthermore, we can find a sequence {vk}k<£N of corresponding eigenfunc- 
tions, that is Vk G -ff^(r2)\{0} and the couple {fJ,k,Vk) *s 0, solution to l)EP2(l for 
any fc G N, so that 

{vk\rA}keN is an orthonormal system of L'^iTA)i 
{llr^jU < — =|rA f *s an orthonormal basis of Hl/^^gf(rA,^)', 



ken 

1 /2 

where we have considered the space H^^j^^f{TA, ^) with the scalar product defined 
in (|^ . 

Beyond the asymptotic behaviour of the eigenvalues, we are interested in the 
asymptotic behaviour of the eigenfunctions, in particular in a kind of exponential 
decay, in terms of the eigenvalues, of the eigenfunctions away from Ta- In the 
next examples, we present some particular cases in which such kind of decay 
holds. 

Example 4.3 Let n = Bn{0, 1) and dn = S^-\ and let Ta = dn and Tj = 0. 



In this case the problem (|EP1|I is a classical Stekloff eigenvalue problem and it 
is well-known that the orthonormal basis of L'^{S^~^) constituted by the traces 
of eigenfunctions, as described in Proposition ^3 coincides with 

(4.3) {fjp : i > and 1 < p < p^} 

where each fjp is a spherical harmonic of degree j , j being a nonnegative integer. 
We have that the function 

(4.4) u,p{x) = \\x\\^f,pix/\\x\\) 



is harmonic in R''^ and solves the eigenvalue problem ljEPl|l with eigenvalue 
A — j. So, the sequence {ujp : j > and 1 < p < Pj} coincides with the 
sequence of eigenfunctions we have described in Proposition 14.11 The integers 
Pj are the dimensions of the spaces of spherical harmonics of degree j and we 
have that, see for instance ^1 page 4], 

f 1 if .7 = 0, 

Pj = \ {2j+N-2)U+N~3y. -f ■ > 1 

( j\{N-2y. li J ^ 

so that 

P,<2(j + 1)^-2, j>0, 

and 

n n 

Ni{n) < J2pj < Xl^O' + 1)^"^ < 2{n+l)^-\ for any n G N. 
j=o j=o 

Furthermore, for any tq, < tq < 1, there exist two positive constants, 
C(ro, N) and a{ro), so that for any Ujp as in H4.4|l it holds 

(4-5) ||wjp||Hi(s«(o.ro)) < C{ro,N) exp{-a{ro)j). 
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Example 4.4 Let il ^ {x e Bn{0,1) : xjy > 0}, and let Ta ^ {x £ 
dBNiO,l) ■■ XN>0} = 5*+-^ andT/ = {x G 5^(0,1) : xn = 0} = B^_i(0,l). 
First of all, we notice that the hypoteses of ProDOsition l4. II and Proposition l4.2l 
are satisfied, so the conclusions of Proposition 14.11 and of Proposition 14.21 hold 
for the eigenvalues and eigenfunctions related to problem (|EP1|I and problem 
(|EP2p with these data, respectively. 

The following exponential decay property can be obtained, as well. 

We have that if w e iJ^(ri)\{0} solves (|EP1|) for a constant A, then, by a 
reflection argument, it follows that there exist j, a nonnegative integer, and /, a 
spherical harmonic function on S^^^ of degree j, so that u{x) = ||a:|p/(2;/||a;|j) 
for any a; € fi and A = j. Thus, if we assume that |l/||L2(rA) = li '^an 
conclude that for any tq, < < 1, 

(4.6) \\u\\m{BN{o,ro)nn) < C(ro, A^) exp(-a(ro)A), 

where the constants C{rQ,N) and a(ro) coincide with the ones obtained in 
Example El 

Again by a reflection argument, we have that if u e -ff^(O)\{0} and a con- 
stant /Li solve (IEP2|I then there exist j, a positive integer, and /, a spherical 
harmonic function on S^~^ of degree j, so that v{x) = \\x\\^ f{x/\\x\\) for any 
X G Q and fx = j. Thus, if we assume as before that ||/||l2(p^-) = 1, we immedi- 
ately infer that for any tq, < tq < 1, 

(4-7) \\v\\h^b M{o,ro)nn) < C{ro, N) exp{-a{ro)n), 

with the same constants C(ro, N) and a{rQ) as before. 

Example 4.5 Let Q = Bn{0, l)\{x S Bn{0, 1) : xn~i > and a;Ar = 0}, and 
let Tj = {x e Bn{0, 1) : Xn-i > and xn = 0}, as in Subsection 12.41 We set 

as in Subsection 12.41 as well. Also the notations concerning Sobolev spaces 
on F^ are the ones introduced in Subsection l2.4l 

The eigenvalue problem ljEPl|) with these data can be rewritten as 

{Aw = in n, 
p^^u ouTa, 
1^ = on (either sides of) F/, 

that is, u e H^{n) solves (|EPr)l if 

/ Vu-Vw— / Xuw, for any ui e if^(ri). 
Jn Jta 

Then, all the conclusions of Proposition l4. II still hold true for the eigenvalue 
problem (|EPl'p . also with the possibility to replace the space L^(Fyi) with the 
space L'^{S^-^). 

The exponential decay of the eigenfunctions is still valid. By separation of 
variables, we have that if m G iJ^(ri)\{0} solves (|EP1'|I with a constant A, then 
there exists a function g G L'^{S^~^) so that u{x) = \\x\\^g{x/\\x\\) for any 
X gVI. Assuming that ||(7||£,2(5n-i') — 1, we obtain that for any rg, < rg < 1, 

(4-8) \\u\\m(B„(ia,ro)nn) < C\{ro, N) exp{-a{ro)X), 
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where Ci(ro, A^) is a positive constant not depending on A and a(ro) coincides 
with the one defined in Example 14.31 



For what concerns the eigenvalue problem l|EP2|) with these data, that is 

Aw = in n, 
(EP2') { ^^^iv ouTa, 

V — onTj, 



we have that v solves (IEP2'(I . in a weak sense, for a constant /i, if u G (^^, P/) 
and 

/ Vw • Vw = / ^vw, for any w G -ffg (SI, P/). 
Jn Jta 

Then, all the results of Proposition 14.21 are still valid for the eigenvalue 
problem (|EP2'|I . and we can again replace the space L'^{Ta) with the space 

If w S Hq{Q., P/)\{0} and /i are a solution to ljEP2'|) . then, by separation of 
variables, we can find a function g e L'^{S^^^) so that v{x) = ||a:||''f/(a;/||a;||) 
for any a; £ $7. If we further suppose ||(7||^2(5jv-i) = 1, we have that for any tq, 
< ro < 1, 

(4-9) ll«llHi(Sjv(o,ro)nn) < C'2(ro,iV)exp(-Q!(ro)Ai), 

where C2(ro, N) is a positive constant not depending on /i and a(ro) is the same 
as before. 



5 Proofs of the main results 

In this section we apply the abstract theorem to the inverse problems described 
in Section 121 and we conclude the proofs of our instability results. 

5.1 Inverse inclusion problem 

We treat the problem related to the Dirichlet-to-Neumann map, the one re- 
lated to the Neumann-to-Dirichlet map and the one related to the experimental 
measurements separately. 



Dirichlet-to-Neumann case 

Proof of Proposition I2. 31 (Dirichlet-to-Neumann case). We apply The- 
orem 13.11 to prove Proposition 12.31 for the operator A. Let us show how the 
abstract theorem is used in this situation. 

We fix xo € X to be equal to Bn{0, 1/2) and we notice that Proposition l2.2l 
implies that X satisfies assumption 0|| of Theorem 13.11 with constants Eq and 
Ci which depend on N, m and [3 only, and with the constant ai = {N — \)/m. 
Furthermore, about X, we remark the following property. We have that every 
-D e X is compactly contained in i37v(0,4/5). 

Concerning assumption[2Il of Theorem l3.1l we observe that for every D £ X , 
the operator A(Z?) is a bounded and linear operator between H^^^(dVt) and its 
dual. We fix if = H^/'^{dVt) and F : X C{H,H') as F{D) = A{D) for 
any D G X. With Fq we denote the Dirichlet-to-Neumann map associated to 
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the problem (|2.8|) - (|2.9I) when D = that is when no inclusion is present, the 
conductor is therefore homogeneous and its conductivity is identically equal to 
1 in f2. Concerning the function 7 our choice is the following. For any ij: G 
i/i/2(afi)\{0}, let 



(5.1) 



7(V') 



11^-112 



Then it remains to choose an orthonormal basis of iJ, {wfe}fceN, so that 7(wfc) 
is finite for any k and (|3.2|l and H3.3() are satisfied. Recalling Example 14.31 in 
particular 1)4.31) . we consider the set 



(5.2) 



3P 



VT+7 



J > and 1 < p < Pj 



with the natural order. This set, by Proposition l4.1l is an orthonormal basis of 
H and it is the one we choose. We also recall that fjp is a spherical harmonic 
of degree j so that \\fjp\\L^(dQ,) = 1, hence 7(/ip/Vl + j) = 1 + J, for any j 
and p. Fixed n G N, #{fc e N : 7(wfc) < n] is clearly bounded from above by 
2(1 + n)^-\ see Example Ol 

We now verify that (|3.3|l is also satisfied. First of all, we remark that for any 
D £ X we have that the operator F{D) — Fq is self-adjoint, in the following 
sense: 

((F(Z?) - Fo)^P, if) = {{FiD) - Fo)^, V), 

for any ip, tp £ H^^^{dfl), where (•, •) is as usual the duality pairing between H' 
and H. In fact, let u be the solution to H2.8|l - H2.9|l and uq be the solution to the 
same boundary value problem with D replaced by the empty set, and let v and 
vq be the solutions to the same boundary value problems with 'ip replaced by (p. 
Then 



{{F{D) - Fo)^, ip)^ / (1 + (a - l)XD)yu ■ Vv 



which is clearly symmetric, thus the self-adjointness follows. 

By self-adjointness, it is enough to show that there exist positive constants 
C2 and a2, depending on N, m, [3 and a only, so that 



(5.3) 



{F{D) ~ ^^o) 



VT+7 



< C2exp(-a2(l+i)), 



for any j and p. In fact, by H5.3|) . we infer that 



([F{D)-F^)^^,^)\ < ||(F(i?)-Fo)J^ 

< C2exp(-a2(l+j)), 



fkq 



n+k 



where we have used the fact that 

self-adjoint, we can reverse the role of j and k and thus obtain 

fjp fkq 



= 1. Since F{D) - Fq is 



((F(i?)-Fo)^^,_= 



< C2exp(-a2(l + A:)) 
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as well, and so 1)3. 3|l immediately follows from 15.3|l . 

Let Ujp{D) be the solution to H2.8|l - H2.9() with the boundary datum ip re- 
placed by fjp, let Ujp be defined as in (|4.4() and let Vjp — Ujp{D) — Ujp. Since 
every D €i X is compactly contained in _Bjv(0,4/5), we can find a constant C3, 
which depends only on N, so that for any D £ X 



1/2 



<C3 / \\Vv,p 

H-^/^(an) \Jo\s„(o,4/5) 



2 



In fact, there exists a constant C3, depending on and 4/5 only, so that for 
every ip G H^/^{dVt) there exists £ with the properties that V'lan — ip, 

EE on Sjv (0,4/5) and 

M\m{n) = IIV'llHi(n\B„(o.4/5)) ^ C'3|lV'||ffi/2(aa)- 
We have that Vjp solves the boundary value problem 



(5.4) 



div((l + (a - l)x£))Vwjp) = -div((a - l)x£)Vujp) in il, 
Wjp = on dfl. 



By the weak formulation of (|5.4() and the fact that D C Bn{0, 4/5), we have 
that for every -tp e H^/^{dn) 

{iF{D)~Fo)f,p,ib) = / Vz;,p-VVi, 
Jo 

thus 



\mD)~Fo)f,p,^)\ < ^j^^.^^-^^\\vv,pry^\j^^ 

^ ^3 (0,4/5) ll^^^jpll'j UWm/Hd^)- 

Again by the weak formulation of H5.4() . we have that 

/o ||Vw,p||2 < max{l, l/a}^l + [a - 1)xd)|! Vi;jp||2 
= max{l, l/a}(— /^(a — l)Vujp • Vvjp) 
< max{l, l/a}|a - 1| (/^ ||V^;,pp)'/' (/^ ||Vu,p|p)'/' . 

From here, using again the fact that D is contained in Bn{0, 4/5), it is easy 
to infer that there exists a constant C4 depending on a only so that 

sl/2 , ^ 111 

\\^v,pf\ <CA\ ||Vu,p"2 



/O\Bjv(0,4/5) I \J Bn{0A/5) 

and so H5.3|) is proved by using 14.5|l . □ 



Neumann-to-Dirichlet case 

Proof of Proposition 12. 31 f Neumann-to-Dirichlet case). For what con- 
cerns the Neumann-to-Dirichlet case, the proposition can be proved as an ap- 
plication of the abstract theorem (see for a similar procedure the proof of the 
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"insulating crack & Neumann-to-Dirichlet case" for the inverse crack problem, 
page r29|) . However, it is also possible to argue in the following way. We observe 
that, for every D N{D) and K{D), the restriction of K{D) to oH^/'^{dn), 
are inverse to each other. Since we have already established the Dirichlet-to- 
Neumann case, there exist ei > and C, depending on N , M, (5 and a only, 
so that for any e, < e < £i, there exists Di and D2 in X so that (|2.11|) and 
(I2.12|) are satisfied. Using the identity 

(5.5) M{D^) - M{D2) = M{D2){k{D2) - MDi))M{Di), 
we infer that 

(5.6) \W{D^)-Af{D2)\\ < \\U{D2)\\fA{D2) - A{D,)\\\\U{D,)l 

where the natural norms have been used. We remark that there exists a positive 
constant C5, depending on N, M, /3 and a only, so that for any D e X we have 

(5.7) \mD)\\ < C5. 

So fTT^ is proved by fTT^ . (gH) and (EH). □ 



Experimental measurements case 

Proof of Proposition 12. 4L For the basic properties of the problem H2.14|l . 
we shall always refer to 22 . Let H = iJ^(ri) x M^. For any D ^ X, and for any 
{u,U) and {w,W) G -ff,"let 

BDi{u,U),iw,W))^ [ (1 + (a - l)xD)Vii • Vu; + V / (u - Ui)iw - Wi). 

•In Jei 

Then, for any 77 G oH^^^^{dV,) and any / G so that X^^^i ^ have 
that there exist, and it is unique up to an additive constant, a couple (u, U) H 
so that we have 

L 

(5.8) BD{{u,U),{w,W))=rj{w\an)+J2^^^^' for any (z«, VF) G iJ. 

1=1 

Furthermore, we have that there exists a constant Cg, depending on A^, a, Z 
and the electrodes only, so that 

(5.9) \\^u\\mn)<Cei\\vlH-^/2^an) + \\I\\)- 
We also remark that if {u,U) solves (|5.8|l . then 

div((l + (a- l)xu)Vw) = infl. 

Therefore, we can rewrite the equation H5.8() on the boundary as follows. We 
denote = f^|aa G QH~^^'^{dfl) and we recall that M{D) is the Neumann-to- 
Dirichlet map associated to the inclusion D. Then, we deduce by straightforward 
computations that (u, U) satisfies H5.8|l if and only if 

(5.10) Ui = + ^ I u, for any / = 1, . . . , L, 
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and the following equation holds in qH ^''^{dfl) 
(5.11) 0.Ei(^(i.)^-^/ 



1=1 ^ ^ 1=1 ^l^'l 

Let JC{D) : ^H~^/^{dVt) ^ oH-^/^{dn) be the operator defined as follows. 
For any oH~^/^{dVt), 

lCiD)<l> = V - ( UiD^ - I N{D^ ] Xer 



We have that JC{D) is a compact linear operator. Since the equation H5.8(l 
admits, up to additive constants, a unique solution, we can infer that (j5.11|l 
is uniquely solvable, therefore the operator Id + JC{D) is invertible, where Id 
denotes the identity operator. Using (|5.9() . this inverse, which we shall denote 
with )C{D), satisfies 

(5.12) \\IC{D)\\ = \\{Id + IC{D))-'\\<Cr, 

where the constant C7 depends on N, a, Z and the electrodes only. 

For any given current pattern /, that is / € M.^ so that ~ we can 

define / = J2f=i (^JeiJ^'^i) oH^^^^{dn). Furthermore, there exists a constant 
Cs, depending on N and the electrodes only, so that, for any / e satisfying 
Ii = 0, we have 

(5.13) miH-^ion) < Cs\\I\\. 

As it is shown in [52], we have that u solves our direct problem (|2.14|) for a 
given current pattern / if and only if (|5.8|l is satisfied with 77 = 0. Therefore, if 
we take / € so that J2f=i = we have that R{D)I = V where, for any 
/ = 1,...,L, 

(5.14) Vi= I J^{D){IC{D)i) + c\ei\, 



where c is a constant which can be computed by imposing the condition that 
^f^i Vi = 0, that is 

5.15) c — —f . 

In order to establish Proposition 12.41 we observe that, by Proposition 12. HI 
we can find constants £1 > and C, which depend on A'', m, (3 and a only, 
so that for any e, < e < £1, there exists D\ and D2 in X so that (|2.11() 
and H2.13(l are satisfied. We show that these inclusions Di and D2 provide us 
with the instability example also in the experimental measurements case. Let 
us evaluate the norm of R{Di) — R{D2). Therefore, we take / £ M.^ so that 
'Y^i^ih — and we evaluate ||(i?(£'i) — R{D2))I\\- We recall that we have 
posed R[Di)[l] = R{D2)[l] = 0. By (|53Hl and lf^l5|l . we have that 

\\{R{D{) - R{D2))I\\ < C4M{D,)(t{D,)i) -U{D2){IC{D2)i)hHan), 
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where Cg depends on N and the electrodes only. Thus, we can find a constant 
Cio, depending on N and Cg only, so that 



\\{R{D,) - R{D,))I\\ < C,o\mD,){IC{D,)-IC{D2))i\\H^,.ian) + 

and, by H5.7|l . H5.12(l and (|5.1^^(l . we also deduce that 

, imo,) - RiD,))I\\ < C,CsC,omD,)-ICiD2)\\\\I\\ + 

^ ' CrCsCio\mD2)~Af{Di)\\\\I\\. 

It remains to evaluate the term ||/C(Di) — /C(I?2)||- We proceed as follows. 
Using an identity analogous to (|5.5|l applied to the operators Id + K, and /C, and 
recalling H5.12|l . we obtain that 

(5.17) \\JC{Di) - IC{D2)\\ <C^\\JC{Di)-JC{D2)\\<Cn\mDi)-Af{D2)\\, 

where Cn depends on N, a, Z and the electrodes only. 

And so the conclusion immediately follows by coupling (|5.16|l with (|5.17|l 
and using (EUSJ- □ 



5.2 Inverse crack problem 

The proof of Proposition 12.51 follows directly from the abstract theorem stated 
in Theorem 13. II We just need to check that all the hypotheses of Theorem 13. II 
are satisfied. Therefore, the proof is divided into two steps, each corresponding 
to one of the hypotheses of Theorem 13. II 

Proof of Proposition 12. 51 - First step. First, let xq <E X he B'j^_^{0, 1/2). 
Then, by Proposition 12. 21 X satisfies assumption Ojl of Theorem 13. II with con- 
stants £0 and Ci depending on N, m and /3 only, and constant ai — [N — l)/m. 
We recall also that a C 3^(0, 4/5) for any a £ X. □ 
For what concerns the second step, we turn our attention to assumption [Izjl 
of Theorem 13.11 Each case, corresponding to operators and A/i, i = 1,2, 
should be treated separately. We limit ourselves to two cases, namely the cases 
corresponding to Afi and A2 , in order to show the main points of the proof, and 
we leave the details concerning the other two cases to the reader. 



Insulating crack & Neumann-to-Dirichlet case 

Proof of Proposition 12.51 - Second step (Insulating crack & Neu- 
mann-to-Dirichlet case). First, we notice that w is a solution to (|2.16(l - (|2.19(l 

if and only if w G H^{n\a), Jg^ u = 0, and 

/ Vu ■ Vw — ri{w\on), ior a,ny w e H^{D,\a). 
Jn\<7 

We observe that for any a G X, Mi{(t) is a bounded and linear operator 
between o-ff~^^^(9r2) and its dual. Hence we take H to be aH-^l'^{d^l) and 
F : X^ C{H, H') to be defined as F{a) = A'i(cr) for any a e X. With Fq we 
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denote in an analogous way the Ncumann-to-Dirichlet map related to (|2.16|) - 
(|2.19|l when cr = 0, that is the Neumann-to-Dirichlet map associated to the 
body where no crack is present. For any 77 e QH^^^^{dQ)\{0}, we define 



(5.18) 



^llH-i/2(Sf2) 



Referring to Proposition l4.1l Example 14.31 and H4.3|l . {wfejfeeN, the orthonor- 
mal basis of H we shall employ, is given by 



(5.19) 



jff. 



J > 1 and 1 < p < I 



with the natural order. We have that j{^/]fjp) — j, for any j and p. Again by our 
remarks in Example 14. 31 we deduce that ^{k G N : 7(wfe) < n} < 2(1 



\N-1 



for any n £N, 

For what concerns H3.3() . we argue in this way. We need a kind of self- 
adjointness of F{a) — Fq for every a £ X. We have that 

{{F{a) ~ Fo)v,cl)) ^ {{F{a) ^ Fo)cl>,ij) 

for any rj, (/) £ oFt^^^'^idfl), where (•,•) is again the duality pairing between 
H' and H. In fact, if u solves (|2.16|l - (|2.19|l . uo solves the same boundary value 
problem with a replaced by the empty set, v and vq solves the same boundary 
value problems with rj replaced by then 



{{F{a) - Fo)v, 



By the self-adjointness of the operator F{a) — Fq, for any a E X, in order 
to prove (|3.3(l we have to show that there exist positive constants C2 and a2, 
which depend on N, m and (3 only, so that, for any j and p. 



(5.20) 



{F{a)-Fo)^f, 



3P 



< C2exp(-a2j). 



We can find a constant C3, depending on N only, so that, for any a G X, 



{F{a)-Fo)y^f, 



3P 



< C'3||«jp||//i(f2\B„(o4/5)); 



where Vjp satisfies 



r Awjp = 



(5.21) 







in f2\iT, 
on 9n, 
on da. 



0, 



with Ujp given by formula (|4.4|l . Since J^^^ Vjp = 0, a Poincare type inequality 
implies that there exists a constant C4, depending on N only, so that, for any 
cr e X, we have 



(i^(a)-i^o)/?/: 



3P 



< C4 



1/2 



O\Biv(0,4/5) 



JPll 
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We can estimate the right hand side of the last equation as follows. We fix a 
cut-off function X so that x G C^(5jv(0, 5/6)), < x < 1, X = 1 on Sjv(0,4/5). 
Without loss of generality, we can assume that for every x G M.^ , || Vx(2:)|| < C5, 
C5 being a constant depending on N only. Let us observe that H5.2HI means that 
for every w € H^{n\a) we have 

/ Wvjp-Ww^- j^^/Vm^p • V(xw). 
Jn\a- Jn\<T 

Then, by taking w — Vjp, we infer that 

Jn\cr JO\cr 

Straightforward computations allow us to prove that there exists a constant Cg, 
depending on TV only, so that 




Then we can conclude using H4.5|l . □ 



Conducting crack & Dirichlet-to-Neumann case 

Proof of Proposition 12. 51 - Second step (Conducting crack & Diri- 
chlet-to-Neumann case) . We begin with a description of the weak formula- 
tion of the boundary value problem H2.17|l - H2.18|l . With iJ^Q„j,j(f2, a) we denote 
the subspace of iJ^(r2) functions which are constant on a. For any c € R, we 
set H^{n,a) as the subset of iJ^(f2) functions which are equal to the constant 
c on a. For any c S M, we have that there exists and it is unique a solution to 
the following boundary value problem 

{Auc = in ^\a, 
Uc = c on a, 
Uc = ip on 951, 

that is a function Uc G H^{Q,, a) so that udon = 4' and that 




Vuc • Vw = 0, for any w € a). 



Given Uc, solution to (|5.22|l . we can define ^\dn e H-^l^{dVt) as follows 

(^|af2,</')= / Vuc-V<^, 
Jn\a 

where (p G H^^^(dfl) and if> is any iJp (fi, a) function so that 'f>\dn — 

We claim that there exists a unique c € M so that (^flgsi, 1) =0, that is 

existence and uniqueness of a solution to (|2.17|) - (|2.18|) . 

We have that u solves (|2.17|l - (|2.18|) if and only if u e Hl^^^^{^,a) so that 

""Ian — Ip a-nd that 




Vw • Vw = 0, for any w G H^in) n Hl^^^^ifl, a). 
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If we take tp to be any HQ{fl,a) function so that ip\dn ~ ipj we have that u 
solves (|TT7|) - (|TTH|) if and only if {t = u - -0 belongs to H^{n) n H^g^^t{n,a) 
and satisfies 



• Vw, for any w € H^{n) n d). 



Standard elliptic theory provides us with existence and uniqueness of such a 
solution. By the property (f^|an, 1) = 0, we can infer that f^|ao G 
can be also defined as 



0\cr 

where ip G H^/^{dVL) and is any _ff^Q„^((r2, cr) function so that = 

Now we can denote with H the space H^^^{dQ), and we can fix 7 as in 
(|5.1I) and the orthonormal basis as the one described in (|5.2|l . The map F : 
X 1-^ C{H,H') is given by F{a) — A2(cr), for any a G X, and Fq denotes the 
Dirichlet-to-Neumann map corresponding to cr = 0. We recall that the operator 
F{a) is self-adjoint for any cr S X, as well as Fq is. 

We proceed to verify in this case. First, there exists a constant Cr, 

depending on N only, so that, for any a Cz X, 



(5.23) 



(Fia) - Fo) 



VT+7 



H-i/2(an) 



1/2 



r!\Siv (0,4/5) 



JPl 



where Vjp — Ujp{a) — J'^'^. , Ujp{a) being the solution to H2.17|l - H2.18|l with 



replaced by j^^. and Ujp being as in H4.4|l . 



Hence, satisfies 



(5.24) 



v„ ^ 



c — 



^lao,l)-0, 



in ri\(T, 
on 9r2, 
on 9(7, 



where c = Ujp{(j)\„. We notice that, if x is the cut-off function previously defined 
in this subsection, then Wjp = [vjp — c + X ^'1+ j ) ^ ^0 (^j Wjplan = — c. 

So, 

/ Vwjp • Vw.jp = (^p^lan, -c) = 0, 

/0\(T 



that is 



/ Vwjp • Vwjp = / Vvjp ■ V(x 



ijp 



from which we easily deduce that 



(5.25) 



1/2 



Sl\o- 



< Cs\\Ujp\\Hi{Bt,{0,5/6), 



where Cs depends on N only. 

So (|3.3|l is obtained by combining H5.23|l . (|5.25|) and (|4.5|l and the self- 
adjointness of the operator F{<t) — Fq. □ 
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5.3 Inverse cavity problem, inverse surface crack problem 
and inverse boundary material loss problem 

As we have already observed, the inverse problem of cavities can be treated 
in a way which is completely analogous to the treatment of the inverse crack 
problem. 

Proof of Proposition 12.61 It can be obtained along the Hnes of the proof of 
Proposition l2.5l with obvious modifications. In particular, the reference point in 
X is given by ao, the orthonormal basis used are those described in Examplc l4.5l 
whereas the reference operator is the one related to the domain il, as in 

Example IHni □ 
Proof of Proposition 12.71 Also the arguments for the proof of Proposi- 
tion l2.7l are simple modifications of what we have used to prove Proposition l2.5l 
clearly making use of the orthonormal basis described in Example 14.41 □ 



5.4 Inverse scattering problem 

The proof of Proposition 12 . 81 is somehow different from the proofs of the analo- 
gous propositions discussed previously. In fact, we can not prove Proposition l2.8l 
as a straightforward application of Theorem l3.ll Nevertheless, the procedure de- 
veloped during the proof of Theorem 13 . 1 1 can be adjusted in such a way to cover 
also the inverse scattering case framework. In the sequel, we limit ourselves to 
the sound-soft case, the sound-hard case can be obtained with minor adjust- 
ments. We shall point out the main differences between the sound-soft and the 
sound-hard case and conclude the proof for the sound-hard case at the end of the 
subsection. The proof of Proposition 12 . 81 for the sound-soft case will be divided 
in two steps. 

Proof of Propositio nI2.8I- F irst step (Sound-soft case). First, we fix 
xo & X to be equal to i?Ar(0, 1) and we observe that assumption 0) of Theo- 
rem lH.ll is satisfied, by Proposition l2.2l with constants Eq and Ci, depending on 
N, m and /3 only, and constant ai — {N — \)/m. □ 
The second step deals with the main difference from the previous cases, 
which is as follows. In the abstract theorem, we have a function F which maps 
elements of a metric space X into elements of C{H,H'), H being a separable 
Hilbert space. Now, fixed a > 0, we define a map F which associates to each 
obstacle D Q X a complex- valued function defined on S'^'^ x S'^~^, namely 
F{D) = As{D){-,-,a) or, respectively, F{D) = Ah{D){-,-,a). In the abstract 
theorem, fixed a suitable Fq G C{H, H'), the operator F{x), x ^ X, was charac- 
terized by the numbers bk,i = {iF{x) — FQ)vk, vi), where /c, ^ G N and {wfcjfegN is 
a suitably chosen orthonormal basis of H . The fundamental properties of such 
a characterization were summarized in assumption of Theorem 13.11 In par- 
ticular, the crucial property was a control on the asymptotic behaviour of the 
coefficients hk^i, which was provided by formulas (|3.2|l and H3.3|) . We shall obtain 
a completely analogous characterization by decomposing the far-field pattern in 
spherical harmonics. 

Proof of Proposition 12.81 - Second step (Sound-soft case). We take 
{ffejfcgN as the orthonormal basis of L?'{S^^^) described in Example 14.31 pre- 
cisely in H4.3|l . with the natural order. Therefore, for each fc G N, Ufe is a (real- 
valued) spherical harmonic function on S^^^ . We set 7(wfe) as the degree of the 
spherical harmonic function Vk- We have that 7(wfe) is an increasing sequence. 
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with respect to k, whose asymptotic behaviour satisfies (|3.2|l with constants 
Ca =2 andp = iV-l. 

The decomposition of tlie far-field pattern in spherical harmonics is given 
by, for any (x,w,a) e S^-^ x S^-^ x (0,oo), 

(5.26) AsiD){x,uj,a) = ^bkjia)vkix)vi{uj), 

k,l 

where the complex- valued coefficients bk,i{a) are given, for any a G (0, oo), by 



(5.27) bk,i{a) = // AsiD){x,uj,a)vk{x)vi{u})dxdLj. 
Furthermore, we use the following characterization 

(5.28) bk.i{a) = / bk{uj,a)vi{uj)du!, 

where the complex-valued coefficients bk{u!,a) are, for any cu E S^^^ and 
any a £ (0, cx)), the Fourier coefficients, with respect to the orthonormal ba- 
sis {wfcjfceN, of the far-field pattern u^(-;a;;a) corresponding to the scattered 
field of the solution to lfT^ - lfOHl - |[?3^ . that is 

(5.29) bk{uj,a)^ / As{D){x,uj,a)vk{x)dx. 

In the next lemma, we establish the asymptotic behaviour of the coefficients 
bk^u which will play the role of the assumption stated in 



Lemma 5.1 Under the previous assumptions and definitions, there exist posi- 
tive constants C2 and a2, depending on N, m, [3 and In only, so that for any 
D G X , for any a £ In and for any (k, ^) € N x N, we have 

(5.30) \bkA")\ < C'2exp(-Q;2max{7(wfc),7(w;)}), 
coefficient bkj as in (|5.27|l . 

Proof. First, we claim that there exist positive constants C3 and a^, depending 
on N, m, /3 and In only, so that for any D E X, for any a E In, for any lu E S^~^ 
and for any fc G N, we have 

(5.31) \bk{uj,a)\ < C3exp(-Q;37(ufc)), 

bk defi ned by fO^. 

By (|5.31|l and H5.28|l . we immediately infer that, for any k, I E N, 

(5.32) \bkAa)\ < |5^~V/'C3exp(-a37(^^fc)), 

|^Af-i| i-,gjj^g (^jy _ l)-dimensional measure of 5^^^. 

Then, we make use of the following reciprocity relation, see for instance |10l 
Theorem 3.13]. For any D E X and any a E (0, 00) we have 



(5.33) As{D){x,uj,a) = As{D){—uj,—x,a), for a.ny x,uj E S 



N~l 
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The reciprocity relation plays the role of self-adjointness for the elliptic operators 
we have considered before and allows us, using H5.32|l . to easily conclude the 
proof of the lemma. Therefore, what remains to be proven is the claim in H5.31|l . 

In order to prove H5.31|l . we begin with a uniform bound on the scattered 
field. We notice that, for any D G X , D C Bn{9/5). With a procedure which is 
analogous to the one used first in |14[ Lemma 2] and later in |2DJ , and using the 
fact that the scattered fields are radiating solutions to the Helmholtz equation, 
we can find a constant C4, depending on N, m, (3 and /jv only, so that, for any 
D G X , any oj e S^~^ and any a e /^r, we have 

(5.34) \u''{x\uj\a)\<Ci\\x\\-'^^-^'^'^, for any a; G M^\Bjv(0, 2), 

where is the scattered field corresponding to the solution to H2.29|I - (I2.3()|I - 
(|2.32|) . We remark that only in the estimate above the difference between the 
cases N = 2 and N = 3 shows up. We refer to 1201 for a detailed discussion about 
uniform estimates of decay at infinity for radiating solutions to the Helmholtz 
equation. 

By Theorem 2.14 in rTUj, we have that since is a radiating solution to 
the Helmholtz equation, with coefficient a > 0, in R.^\Bn{9/5), then, for any 
X e R^\BAr(0,2), 

(5.35) u%x;uj;a)=Y,bk{uJ,a)H^f^jV^\\x\\)vk{x/\\x\\), 

fc 

where bk are complex-valued coefficients given by 

i^M) ».<.,(v^.)^/^___..(.*;.;aK(i)cl*, to, any .> 2. 

where, for any integer n > 0, Hn'^ denotes the Hankel function of first kind and 
order n. 

Theorem 2.15 in 10 provides us with the necessary link between coefficients 
bk and bk- In fact, it holds that 

(5.37) hiu;, a) = (V2)(^-^)/^a-(^-i)/4(-i)^''+(^-i)/^6,(^, a). 

We combine (|5.37ll with H5.36|l and, by using H5.34|l . we obtain that there 
exists a constant C5, depending on C4 only, so that, for any fc e N, any u £ S^~^ 
and any a In, 

(5.38) \kiLu,a)\ < C5a-(^-i)/V-(^-i)/2|ijW^^(yHr)r\ for any r > 2. 

We choose r as follows. For any a, a > 1, we take r = 2, whereas for any a, 
< a < 1, we pick r — 2/ -/a. With this choice we infer that for any a e /at, we 
have that 

(5.39) \k{L^,a)\<Ce\H^ljP)\-' 

where 2 < r < 2max{l,a} and the constant Cg again depends on C4 only. 
Then we can establish H5.31|l by using the well-known asymptotic behaviour of 
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the Hankcl functions. In fact, there exists a constant C7, depending on N and 
a only, so that for any f, 2 < f < 2inax{l,a}, 

(^•40) ^ { S (f )" in - !)-("-) LTny 2. 

By a straightforward computation, H5.39|l together with (|5.40() imphes the va- 
lidity of (ESU). □ 
Now we have what is needed to prove ProDOsition l2.8l in the sound-soft case. 
Let us just notice that, for any D € X and any a € (0, 00), 

IIA(i?)(-,-,a)||L2(s«_ixs«-i) = I ^|6fe,i(a)|2 

Then, with the same procedure used to prove Lemma l3.2l and keeping in mind 
the fact that we repeat the procedure j times, one for each a;, i = 1, . . . , j, 
we can find a constant Cg, depending on N, and the constants C2 and 012 of 
Lemma [5.11 onlv. so that, for any S, < S < 1/e, there exists a subset Yg of 
X with at most exp{jCs{—^ogS)'^^^'^) elements so that for any D G X there 
exists D G Ys satisfying 

sup AP))(-,-,a)||L2(s«-ixs"-i) < 

Then the conclusion of the proof of Proposition 12.81 in the sound-soft case is 
immediate. □ 
Proof of Proposition 12.81 (Sound-hard case). We conclude this subsec- 
tion sketching the proof for the sound-hard case. 

First, for any a, < a < 1, and /3 > we define Y(i,a)^(i/2)(<S'^^^(0, 1)) as 
we have defined F„^(i/2) {S^^^{0, 1)), m being an integer, with the only obvious 
modification of replacing the C™ norm with the C^'" norm. Furthermore, we 
observe that y(i,Q)/3(i/2)(*S'^^^(0, 1)) satisfies assumption[3l of Theorem l3.1l with 
constants Eq and Ci, depending on N, a and (3 only, and constant ai — (N — 
l)/(l + a). 

The difference between the sound-soft case and the sound-hard case relies 
in the estimate contained in H5.34|l . With arguments which are analogous to 
the ones used for the sound-soft obstacles, estimate (|5.34|) can be proved for 
sound-hard obstacles belonging to i^m/3(i/2)(5'^^^(0, 1)), for any integer m > 2, 
with a constant C4 depending on N, m, [3 and /jv only, and for sound-hard 
obstacles belonging to y(i,Q)^(i/2)(5'^~^(0, 1)), with < a < 1, with a constant 
C4 depending on N ^ a, (3 and Im only. 

Since the other part of the proof does not depend on the type of boundary 
conditions used, the result follows, for sound-hard obstacles, for any m > 2. 
Some modifications are needed to treat the case when m = 1. We have that, 
for any < a < 1, Y(i,„)^(i/2) (5^-1(0, 1)) C Y^p^^,2){S'''\Q,l)). We apply 
the procedure described before to the set Yi^i^a.)i3(i/2){S'^~'^{Q, 1)) and, recalling 
Remark 13.31 we infer that for any 7 > there exists a constant £1 > 0, that 
depends on iV, j, /3, /at, a and 7 only, so that for any e, < e < £1, we can find 
D3 and Di, both belonging to y(i,Q)^(i/2)('5'^~^(0, 1)), satisfying ^(£'3, D4) > e 
and 

sup ||(A(-D3)- Ap4))(-,-,a)llL2(s«-ixS"-i) < 2exp(-£""i), 

a£{ai,...,aj} 
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where ai = 



N-l 



. We can choose, from the very beginning, a and 7 



(l+a)(2Ar-l+7) 



in such a way that a and 7 depend on N only and (1 + a){2N — 1 + 7) = 2N, 
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